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ABSTRACT 



We present the results of a near-infrared imaging study of high redshift (z ~ 3) quasars using the ESO-VLT. Our targets were selected 
to have luminosities among the highest known (absolute magnitude Mb ;SS —28). We searched for resolved structures underlying 
the bright point-source nuclei by comparing the QSO images with stars located in the same fields. Two QSOs (HE 2348—1444 at 
z — 2.904 and HE 2355—5457 at z = 2.933) are clearly resolved in K s , and with somewhat lower significance also in H\ one object 
is resolved only in K s . At these redshifts, H and K s correspond almost exactlly to rest-frame B and V, respectively, with virtually no 
if -correction. We also report briefly the non-detection of some additional QSOs. The detected host galaxies are extremely luminous 
with Mv ~ —25. Their rest-frame B — V colours, however, are close to zero in the Vega system, indicating substantial contributions 
from young stars and a stellar mass-to-light ratio below 1 (in solar units). Tentatively converting My and B — V into rough estimates 
of stellar masses, we obtain values of M* in the range of several 10 11 Mq, placing them within the high-mass range of recent 
high-redshift galaxy surveys. We present optical spectra and use C IV line width measurements to predict virial black hole masses, 
obtaining typical values of Mbh ~5x 10 9 Mq. With respect to the known correlation between host galaxy luminosity Lv'.host and 
Mbh, our measurements reach to higher luminosities and redshifts than previous studies, but are completely consistent with them. 
Comparing our objects with the local (z ~ 0) Mbh-Afbuigc relation and taking also the low stellar mass-to-light ratios into account, 
we find tentative evidence for an excess in the Mbh/A^buigc mass ratio at z ~ 3. 

Key words, galaxies: active - galaxies: evolution - galaxies: high redshift - quasars: general 



1. Introduction 

Our understanding of quasar host galaxies has made rapid and 
significant progress over the past decade. It is now clear that lu- 
minous quasars (Mb Si —24) at low redshifts are mostly hosted 
by m assive elliptical galaxies dMcLure et alJll999tlDunlop et alj 
2003), irrespective of r adio power. Nuclear and host bulge lumi- 
nosities are correlated (M cLeod & Riekel [T995). which is natu- 
rally explained by the tight relation betwe en the masses of su- 
permassive black holes and bulges (e.g., lHaring & Rixl I2004 
and references therein), in particular if quasar black holes typ- 
icall y accrete and radiate at a fixed fraction of the Eddington 
rate (Mc Leod et aTlll999t) . 

High redshift quasar hosts have also received increased at- 
tention lately, as they provide one avenue to study the assem- 
bly and evolution of massive galaxies, in particular in relation to 
the growth of their central black holes. Yet, quantitative knowl- 
edge of QSO host properties at redshifts well beyond z <; 1 
is still very limited, mainly due to the fact that optical imag- 
ing only gives access to the rest-frame ultraviolet spectral range 
where the QSO nuclei outshine any host galaxy. The improve- 
ments in near-infrared instrumentation have ope ned also this 
window, first with the Hub ble Space Telescope (Kukul aet alj 
1200 It iRidewav et all l200ll). fol lowed by some ground-base d 
experiments dFalomo et alj|2004 120051: |Kuhlbrodt etaT] 120051) . 
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Very useful insights have also b een obtained from HST imaging 
of gravitationally lensed QSOs dPeng et alJl2006bl) . 

These studies were all limited to single band observations, 
typically sampling the rest-frame optical wavelength range. 
Largely focussed on solving the challenging problem of resolv- 
ing the QSO images and deblending the nuclear and underlying 
extended galaxy components, their main result is that QSOs are 
hosted by fairly luminous gala xies already at high z. More quan- 
titatively, [Kukulaelal. (2001) found essentially no evolution be- 
tween z ~ 2 and z ~ in typical host luminosities, for radio- 
quiet QSOs at given absolute nuclear magnitude of Mb ~ —24. 
This trend has been confirmed by mo re recent studies based on 
larger samples (e.g. Peng et al. 2006a]b|). 

Unfortunately, additional information such as morphologi- 
cal properties or colours are even harder to come by than single 
band luminosities. Especially colours, or more generally speak- 
ing, spectral data are very scarce already at much lower red- 
shifts, and significant improvemen ts have only been obtained 
in the las t few years, at z ^ 0.2 dJahnke et alj l2004ah and at 
z ~ 0.7 (Sanchez et al. 2004). These studies have demonstrated 
that there is strong evidence for the presence of of blue, young 
stars in low-redshift QSO host galaxies, including those that ap- 
pear to be quiescent ellipticals. 

Clearly, one would expect recent or ongoing star formation 
to be even more prevalent in galaxies at high redshifts. Indeed we 
found substantial UV emission from the hosts of mod erate lumi- 
nosity QSOs at z ~ 2.2 using HST/ ACS imaging dJahnke et al.l 
2004b). But to our knowledge, no measurements of rest-frame 
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optical colours have yet been published for high-redshift QSO 
host galaxies. 

In the present paper we report on the results of a pilot 
study, aiming at measuring host galaxy luminosities in two dif- 
ferent near-infrared bands using the ESO-VLT. Our targets are 
radio-quiet, very high luminosity QSOs at z ~ 2 and z ~ 3. 
We briefly describe our strategy to decompose the QSO im- 
ages into their nuclear and host contributions, and we discuss 
the results in the context of recently obtained scaling relations. 
Throughout this paper we assume a flat cosmology with Hq = 
70km s" 1 Mpc" 1 , Vt M = 0.3 and A = 0.7. All magnitudes 
are expressed in the Vega system. 

2. Observations and data reduction 

2.1. Target selection 

Our aim was to observe quasars with the highest possible lumi- 
nosities, which we expect to reside in the most massive bulge- 
dominated gal axies. Using the Ham burg/ESO Survey for bright 
QSOs (HES; IWisotzki et aLll2000i Wisotzki et al., in prep.) as 
input catalogue matched this aim. Our QSOs were selected to 
have redshifts 2.6 < z < 3.0 so that the rest-frame B and 
V bands were approximately reproduced by the standard near- 
infrared filters H and K s , respectively. We also added a few 
QSOs with 1.8 < z < 2, where J and H correspond to rest- 
frame B and V, respectively. Figure[T]illustrates the good match 
between rest-frame and observed-frame filter bands. We verified 
by visual inspection in the Digitized Sky Survey that there was 
at least one star of comparable brightness within ~ 1 arcmin 
from each QSO, which later would be usable for calibrating the 
point-spread function. Table Q] lists the small sample of objects 
for which we obtained useful data. 

2.2. Observations 

The sample was observed with the ESO-VLT in service mode 
(thus spreading over several nights) between May and July 2002. 
The instrument used was the Infrared Spectrograph and Array 
Camera (ISAAC), mounted on the Nasmyth B focus of UT1 
(Antu) of the VLT. Because of technical problems with the short 
waveband ('Hawaii') arm, only the long waveband arm was 
available. This was equipped with a SBRC Aladdin 1024 x 1024 
pixel InSb array providing a pixel scale of 0'.'148 pixel -1 for a 
total field of view of 151" x 151". The data were obtained in ob- 
serving blocks consisting of sequences of 12 exposures in J, 15 
in H and 30 in the K s band. In order to account for the rapidly 
varying NIR night sky, the telescope was moved between succes- 
sive exposures following a multi-point dither pattern along a se- 
quence of predefined quasi-random grid points. Exposure times 
per dither pointing were 150, 120 and 60 s in J, H and K s , re- 
spectively, each split into several short integrations. The final net 
integration times per target and filter band are listed in Table Q] 
along with the effective seeing measured in the coadded images. 
All observations were obtained under photometric conditions. 

2.3. Data reduction 

The images were reduced using ESO-MIDAS and the eclipse 
software package also provided by ESO, in particular the 
eclipse-based ISAAC reduction pipeline recipe lw j itter. The 
basic steps involved (i) flat fielding using twilight flats, (ii) esti- 
mate of a sky frame by median filtering of the individual frames 
in the stack, after rescaling to account for temporal sky level 




Fig. 1. Position of the corresponding filter profiles at z = 2 and 
z = 3. In both panels, the spectral template represents a single 
stellar population model of a 1 Gyr old galaxy. 



variations, (iii) sky subtraction, (iv) determination of the angular 
offsets between frames, using an auto-correlation technique, (v) 
final image registration and coaddition for each stack of images. 

As all observations in the NIR, our data are totally domi- 
nated by the sky background. It turned out that the background 
removal by the reduction pipeline was not quite good enough 
for a detailed analysis of very low surface brightness features. 
We therefore determined a local background correction for each 
object separately by creating a curve of growth. This gener- 
ally worked satisfactorily, although the data sometimes suffered 
from spatial background features varying with time, which we 
could not correct for. The main effect of such variations was 
an increase of the statistical uncertainty of the local background 
correction. We quantified this uncertainty by extracting several 
small control images of blank sky in arbitrary positions, which 
were then treated in the same way. 

2.4. Photometry and absolute magnitudes 

Flux calibration was bas ed on observations of standard stars 
from Persson et al] dl998l) . Atmospheric extinction coefficients 
were taken from the ESO web pages. The derived zero points 
from various stars turned out to be highly consistent, with a me- 
dian internal error of 0.003 mag. 

Conversion of apparent into absolute magnitudes requires 
the knowledge of a X-correction (including the band pass 
stretching). By design of our experiment, the spectral corrections 
are very smal l (cf. Fig.[T|). We follo wed in principle the approach 
described by Jahn ke et al.l (l2004al) . For the K-correction of the 
nucleus we used an average quasar SED, for the host galaxies 
we calculated K-corrections from single stellar population mod- 
els with stellar ages ranging between 0.1-1 Gyr. The resulting 
absolute magnitudes, however, are almost independent of the as- 
sumed spectral energy distributions. 

3. Analysis 

3. 1. Point spread function 

The detection of faint extended fuzz underlying a bright point 
source is notoriously difficult and requires foremost an accurate 
knowledge of the effective point spread function (PSF) of the 
optical system. All quasars in the sample were selected to have 
at least one suitable PSF star in their vicinities. Unfortunately, a 
detailed analysis of stellar profiles within the full ISAAC field 
of view revealed significant spatial variations of the PSF, rem- 
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source 


RA 


Dec 
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B 


exp. time [min] 


seeing [arcsec] 














J 


H 




J 


/-/ 




HE 2329- 


-4115 


23 31 50.5 


-40 58 53 


1.850 


18.7 


60 


30 




0.75 


0.70 




HE 2213- 


-3722 


22 16 30.8 


-37 07 32 


1.920 


17.6 


30 


60 




0.71 


0.73 




HE2149- 


-0436 


21 51 50.1 


-04 22 19 


2.643 


17.3 




60 


60 




0.58 


0.71 


HE 2355- 


-5457 


23 58 33.5 


-54 40 42 


2.904 


17.8 




90 


90 




0.73 


0.70 


HE 2348- 


-1444 


23 51 29.8 


-14 27 57 


2.933 


17.3 




60 


60 




0.65 


0.65 



iniscent of optical aberrations. Even worse, the pattern clearly 
varied also in time, possibly due to focus variations, and not two 
images shared the same pattern. While in principle it is possi- 
ble to mod el such distortions usi ng several stars in the field of 
view (e.g., iKuhlbrodt et al.ll2004l) . the number of stars at these 
high Galactic latitudes was mostly too small. Some of our im- 
ages thus became essentially useless, with the PSF stars display- 
ing a broader radial profile than the QSOs. Such data had to be 
discarded entirely, and they are not discussed any further in this 
paper. 

In the end, we used only exposures where the PSF appeared 
to be reasonably stable within the central regions. In our best 
cases we had two stars close to the quasar that were bright 
enough to be useful as PSF calibrators; here we used either 
an average PSF for each given image, or each star individually 
when explicitly quantifying the effects of PSF uncertainty. In 
particular we adopted the difference image of the two stars as 
a rough guess of the additional 'noise' contributed by the PSF 
uncertainty. 

3.2. PSF subtraction 

PSF subtraction is a well-established technique for the analysis 
of QSO host galaxy images. It has the big advantage that a min- 
imum of assumptions about the nature of the host is entering 
the analysis. Here we employed it with the main purpose of per- 
forming a conservative test for the 'null hypothesis', namely that 
no resolved host galaxy is detectable underneath the QSO point 
source. 

Once a PSF was established, it had to be shifted to the lo- 
cation of the QSO, scaled, and subtracted. To minimize shift- 
ing residuals we used a subpixel grid for the centroid position 
and found the best location from \ 2 minimisation. As a scal- 
ing criterion we demanded that the summed residual flux in a 
small aperture of 2 pixels radius around the centroid vanished 
after PSF subtraction. Unresolved sources can thus be robustly 
identified, as any residual flux after PSF subtraction should be 
consistent with PSF mismatch errors. On the other hand, it is 
obviously a very conservative criterion for resolved sources im- 
plying a significant oversubtraction in case a host galaxy should 
be present. However, a significantly detected extended source 
should still yield systematically positive residuals at larger radii. 
Only those detections were taken as significant that show pos- 
itive residuals for most of the bright stars that are located in 
the vicinity around the QSO and contain at least 5 % of the 
QSO flux after oversubtraction correction. In such cases it is 
even possible to approximately account for the oversubtraction 
by applying a corr ection factor derived from simulations (see 
Jahnke_eiaL|20P4bfor a more extensive discussion of this proce- 
dure). Briefly, we performed a series of such simulations where 



artificial host galaxies were added to one of the PSF stars. After 
adding appropriate noise we took a second star in each field as 
PSF calibrator which was shifted, scaled, and subtracted. This 
provided us with a rough estimate of the systematic PSF over- 
subtraction. The model host galaxies in this exercise were as- 
sumed to be smooth and circularly symmetric, with a surface 
brightness distribution described by a de Vaucouleurs profile of 
effective radius of 3, 5 and 10 kpc, respectively. This is slightly 
more compact than present-day high-mass elliptical galaxies, but 
plausible given the high redshifts invo lved and given earlier re- 
sults (see also Kuhlbrodt et al. 2005). While typical values of 
the residual flux are of the order of a few percent we determined 
oversubtraction correction factors of 0.5-1 mag for the K s band 
and somewhat higher values of 1-1 .5 mag in case of the H band. 
We obtained similar correction factors if an exponential disk is 
assumed instead of a de Vaucouleurs profile as galaxy model. We 
find the apparent host galaxy magnitudes to be fainter by 0.2- 
0.5 mag and stronger residuals at the same time. After the over- 
subtraction correction was applied, all nuclear and host galaxy 
magnitudes agree within 0.2 mag with the values determined by 
our two-dimensional modelling technique. A similar procedure 
as described above was applied for the non-detections to deter- 
mine upper limits of the host luminosity. Here we created series 
of quasars images by adding a synthetic galaxy component of 
3, 5 and 10 kpc, respectively, to the given PSF. We varied the 
nuclear/host flux ratio until we clearly resolved the galaxy in 
the radial surface brighness profile. We adopted as upper limits 
when the quasar surface brightness deviated from the PSF pro- 
file by 3(7 based on shot noise error bars determined in a ring 
excluding the central pixels. 

3.3. Two-dimensional modelling 

A more advanced deblending treatment is based on the explicit 
assumption that the surface brightness profile of an observed 
host galaxy should roughly follow a simple analytical prescrip- 
tion. It is then possible to fit a model of the two-dimensional 
light distribution to the data, consisting of a superposition of a 
scaled PSF plus an analytical galaxy mo del convolved with the 
PSF. We used our own code PAMDAI (IKuhlbrodt et al.l [2004) 
which in t he past was successf ully applied to low-redshift QSO 
hosts (e.g. JJahnke et al.l2004al) . The code determines the best-fit 
parameters by % 2 minimisation employing a modified downhill 
simplex algorithm. 

We soon realised that the signal of any possible host galaxy 
was too weak in our data to permit simultaneous constraining of 
the shape and the flux level of the galaxy model. We therefore 
made the following additional assumptions: The host galaxy was 
assumed to follow a de Vaucouleurs law with zero ellipticity and 
with a preset effective radius r c s ■ When assuming an exponen- 
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Fig. 2. Residuals after subtracting one PSF star from another star 
for HE 2355-5457 (left) and HE 2348-1444 (right). We used 
the two stars marked as SI and S2 in Fig. [3] Cuts are chosen the 
same as for Fig. [3] 



tial law instead, we find the host magnitudes to be fainter by 0.3- 
0.4 mag. Thus, only the fluxes of the nucleus and galaxy had to 
be determined by the fit. With these provisions the fits converged 
well for those objects where the PSF subtraction indicated a re- 
solved source. It turned out that the \ 2 values resulting from 
assuming different fixed values of r e g came out very similar as 
long as realistic r e s values were chosen. It is thus not possible to 
significantly constrain the scale length of the host galaxy from 
our data, which is perhaps not too surprising as our angular res- 
olution of 0'.'7 corresponds to ~ 5 kpc at z = 3. However, the 
estimated host galaxy fluxes do not depend se nsitively on the 
corre ct scale length. This fact is known for long ([Abraham et al.l 
fl992h . and we confirm it to be valid for our present data, also 
through extensive simulations. 

We computed a two-dimensional model only in those cases 
where the conservative PSF subtraction revealed positive resid- 
uals. We then adopted the best-fit PSF scaling factor to obtain a 
best-guess PSF-subtracted image of the host galaxy. 



somewhat lower significance) for the H band. However, with the 
galaxy contributing only ~ 5 % of the total uncorrected K s band 
flux, the host is still much fainter than the central nucleus. In the 
H band this fraction is even slightly less than 3 %. After ap- 
plying a statistical correction for oversubtraction (see Sect. 13.21 ) 
these values increase to ~ 10 % and ~ 8 % for the K s and H 
band, respectively. We then computed two-dimensional model 
fits (with fixed scale lengths as discussed above), and obtained 
integrated K s and H band fluxes for the host galaxy that are 
in good agreement with the corrected estimates based on PSF- 
subtraction. This makes HE 2355-5457 the QSO with the high- 
est significance detection in the whole sample. The results are 
documented in the upper panels of Fig. [3] 



HE 2348-1444: The situation with this QSO is similar to 
the previous object: We have two bright stars close to quasar 
which show very similar surface brightness profiles. We again 
find positive PSF-subtraction residuals in both bands which 
were confirmed by the same test procedure as described for HE 
2355—5457. However, the fraction of an extended flux compo- 
nent is even lower, only 1-2 %, or 6 % and 5 % in K s and H, 
respectively, after correction for oversubtraction. Notice that HE 
2348 — 1444 is the brightest quasar in our sample. In this object 
the residual flux is strongly asymmetric (see Fig. [3}. This feature 
becomes particularly evident after subtracting the azimuthally 
symmetric two-dimensional model fit. At the same time, the 
lacking symmetry makes the extended nature of the object un- 
conspicuous in the radial profiles (it is in fact not visible in the H 
band profile). Nevertheless, we believe that since the feature ap- 
pears as highly consistent in both bands, its reality can be safely 
assumed. We conclude that also in this QSO we have detected 
the host galaxy, and that there is evidence for a highly disturbed 
morphological structure, most likely due to gravitational inter- 
action or merging. 



4. Results 

We succeeded to resolve the host galaxies of three QSOs from 
the higher redshift subsample, two of them in both H and K s 
bands, and one only in K s . We failed to detect the hosts of two 
QSOs from the lower redshift subsample. We now briefly discuss 
the sources individually. A summary is given in Table |2j images 
of the detected QSOs are shown in Fig. [3] 

HE 2355—5457: This field has two nearby stars of compara- 
ble brightness as the QSO located within less than 45" distance 
from the QSO. There is no evidence for significant PSF vari- 
ations over the central regions of the field. Subtracting a com- 
posite scaled PSF from the QSO left a strong excess of posi- 
tive pixel values outside of the very central zone, in both H and 
K s bands. We tested for the reality of this effect by subtract- 
ing one PSF star separately from either QSO or the other star. 
The difference of QSO— PSF showed consistently a 10 x higher 
residual flux than the difference of PSF1— PSF2. Visual inspec- 
tion of the PSF-subtracted QSO image showed that this excess 
in the quasar profile was not caused by a close companion, or 
by a strongly asymmetric feature, but that it was due to more 
or less azimuthally symmetric low surface brightness emission 
centered on the QSO. We conclude that our conservative null hy- 
pothesis of the QSO having no detectable host galaxy has to be 
rejected, with high significance for the K s band, but also (with 



HE 21 49—0436: The field contains several (mostly faint) stars 
showing that PSF variations for the whole field of view are rather 
large. For both bands we find positive as well as negative resid- 
uals after subtraction of individual stars extracted within 30" 
around the QSO. In case of the K a band, the number of suf- 
ficiently bright field stars was acceptable to create an analyti- 
cal rnod£lofflie_^atmlPSF variations (similar to the procedure 
bv lKuhlbrodt et al.ll2004b . Using the best-fit analytical PSF pre- 
scription at the QSO position, we detected an extended source 
after PSF subtraction. The amount of uncorrected flux is 6 % 
and increases to 13 % after oversubtraction correction. The im- 
age quality of the H band data was unfortunately inferior to the 
K s data, and we did not succeed in resolving the host galaxy. In 
the following we only quote a 3a upper limit of the host galaxy 
luminosity, propagating into corresponding lower limits of the 
rest-frame B — V colour. 



HE 2329-41 1 5 and HE 221 3-3722: While the J band data 
of these z ~ 2 QSOs turned out to be useless, we found both 
QSOs to be marginally resolved in the H band, but below our 
formal significance criterion. We therefore only quote 3<r upper 
limits in Table [2] 
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Fig. 3. Host galaxy detection of HE 2355-5457 (upper panel), HE 2348-1444 (middle panel) and HE 2149-0436 (bottom panel). 
The left-hand images have a scale of 1' x 1' and show the QSO with PSF calibrator stars; for each QSO, H band in the top row, 
K s band below. All other images are 8" x 8" postage stamps zoomed on the QSO. The three middle panels show the residual host 
galaxy light after PSF subtraction; middle left: based on conservative PSF subtraction; center: based on the best-fit PSF scaling; 
middle right: residuals left after subtracting both the best-fit PSF and an azimuthally symmetric model host galaxy. In the right hand 
panels we show the radial surface brightness profiles of the QSO and corresponding best-fit PSF. 
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Table 2. Integrated apparend magnitudes and absolute magnitude estimation for the observed QSOs. The results are quoted for an 
fixed effective radius of 5 kpc. 



Object 




^^■host 




K a host 


K s nuc 


M nuc 


M host 


M v uc 


M host 


HE 2149 


-0436 


>18.0 


16.3 


18.0 


15.9 


-27.4 


> -25.6 


-27.4 


-25.4 


HE 2355 


-5457 


18.0 


15.8 


17.8 


15.5 


-28.2 


-25.8 


-28.2 


-26.0 


HE 2348 


-1444 


17.8 


14.8 


17.4 


14.2 


-29.3 


-26.4 


-29.5 


-26.3 


HE 2213 


-3722 


>18.5 


16.6 










-26.7 


> -24.8 


HE 2329 


-4115 


>17.8 


16.9 










-26.3 


> -25.4 



5. Discussion 

5. 1. Host galaxy luminosities 

The fluxes of the detected host galaxies were determined by inte- 
grating over the PSF-subtracted images, where the PSF scaling 
factor was given by the best-fit analytic model. Since we did 
not allow the half-light radii to vary during the fit, the resulting 
fluxes depend somewhat on the assumed size of the galaxy. We 
found that our adopted range for r 1 / 2 from 3 kpc to 10 kpc cor- 
responds to a systematic uncertainty of ±0.2 mag, regardless of 
the used filter. In those cases where only upper limits are quoted, 
the conservatively PSF-subtracted images were used. 

After applying distance modulus and K corrections as dis- 
cussed above, we obtained absolute magnitudes in the rest-frame 
V and B bands. We reiterate that our measurements essentially 
sample the objects at these rest-frame wavelengths, and that K 
corrections are almost negligible. The resulting absolute magni- 
tudes are listed in Table[2] Notice that the quasars in our sample 
are so powerful that despite the high nuclear-host ratio of <; 10, 
the hosts are still extremely luminous galaxies. 

Figure [4] compares our measurements to other recent stud- 
ies of high-redshift QSO host galaxies. We restrict this compar- 
ison to near-infrared studies where My of the host galaxy can 
be estimated without major spectral extrapolation. While evi- 
dently there are very few detections for high-luminosity QSOs, 
our data line up well with the published observations. As a whole 
Fig. [4] suggest that QSO nuclear and host galaxy luminosities are 
correlated also at high redshifts, very similar to the finding by 
iFlovd et all d2004l) at low z. 

QSOs radiating at a fixed fraction of the E ddington limit 
are di stributed along a diagonal in this diagram dMcLeod et alj 

1999) , assuming that bl ack hole masses are p roportional to 
host bulge luminosities (Magor rian et al.l Tl998l) . Under these 
premises, our objects appear to have similar Eddington ratios as 
less luminous QSOs at redshifts z 2 (and in fact also similar 
to low-z QSOs). This very qualitative reasoning can be quan- 
tified if mass estimates are available, for both black holes and 
stellar bulges. In the following we attempt to estimate both these 
quantities for our z ~ 3 QSOs. 

5.2. Black hole masses 

Several recent investigations have demonstrated that one can es- 
timate the masses of black holes in active g alactic nuclei from 
single-epoch spectra of broad emission lines dVesterg aard 2002; 
IVestergaard & Petersonll2006t iMcLure & Jarvisll2002l) . This ap- 
proach assumes the bro ad-line region (BLR) to be in approxi- 
mate virial equilibrium ( Peterson & Wandell2000h a nd adopts an 
empir ical luminosity-BLR size scaling relationship dKaspi et"aT] 

2000) . The resulting object-to-object scatter is substantial, but 



-28 



-26 



-24 



-22 



-20 



3) 



mo o o * 
+o 



o o 



EB Ridgway et al. 2002 - 
* Falomo et al. 2004 
O Peng et al. 2006 
+ Kukula et al. 2001 _ 
^ Kuhlbrodt et al. 2005 
■ Current sample 
_l l l l I l l l 



-20 



-22 



-24 



-26 



-28 



-30 



Fig. 4. Relation between nuclear and host galaxy luminosities 
for high-redshift (z > 1.4) QSOs. Our data are represented by 
the filled squares. For comparison we plott several results from 
other recent studies: (jRidgwav et al.|[200l[: [Kukula et al.l [ 200 1 : 
iFalomo et alJ 12004 iKuhlbrodt et ail 120051: iPeng et alJ 12006a). 
all numbers converted to our set of cosmological parameters. 
Objects radiating at constant fraction of their Eddington lumi- 
nosity would lie on a diagonal line in this diagram. 



there is no systema tic bias |Collin et al.ll2006l). W e adopt the 
prescription used by Vestergaard & Peterson! d2006l) to compute 
black hole masses directly from observables, 



M bh = 5.4 x 10 b 



! 


' oi(Civ) " 


2 r 




1000km/s_ 





AL A (1350A) 



10 44 erg/s 



0.53 ' 



(1) 



where L\(1350 A) is the monochromatic luminosity at 1350 A, 
and ct;(Civ) is the dispersion (second central moment) of the 
C IV A1550 emission line. 

From the follow-up spectroscopy within the Hamburg/ESO 
Survey we have slit spectra available of all three QSOs where we 
resolved the host galaxies. These spectra are presented in Fig. 
[5] We measured the C IV line dispersion in the following way: 
We first subtracted a local continuum, estimated as a straight 
line fitted to the line-free windows at 1440 A < A < 1480 A 
and 1660 A < A < 1700 A; a slight extrapolation of this pro- 
vided also the continuum flux at 1350 A. We then determined 
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the second central moment a\ for each line. Two of our objects 
(HE 2348-1444 and HE 2355-5457) show evidence for associ- 
ated absorption in the blue wing of C IV which we had to correct 
for. In HE 2355-5457 it was straightforward to interpolate over 
the absorption line, whereas in HE 2348-1444 some guesswork 
was required, given the low spectral resolution. Fortunately, a\ 
is not very sensitive to such manipulations (unlike the FWHM), 
and the difference between the best-guess corrected and an even 
completely uncorrected o\ was only ~ 10 %, which would trans- 
late into a 20 % change in Mbh- As some degree of correction 
was undoubtedly needed, we argue that the additional error in- 
troduced this way is smaller. The results are presented in Table 

El 

As expected given the high luminosities, the black holes 
are very massive with Mbh U P to a few times 10 9 M Q . Yet, 
these masses are w ell within the upp er envelopes of the distri- 
bution obtained by Veste rgaardl d2004l) for high-luminosity high- 
redshift QSOs. In Fig. [6] we investigate where our objects fall 
within the Mbh — Ly relationship, compared to other high red- 
shift QSOs with detected host galaxies and available black hole 
mass estimates. (In this context we use the term Ly to avoid 
confusion between masses and absolute magnitudes; neverthe- 
less we always express luminosities in terms of absolute mag- 
nitudes My.) Our three objects are perfectly consistent with an 
approximately linear Mbh ~ log Ly relation. 

A revised prescription of the Mbh ~ Ly rel ation for low- 
redshift bulges and ellipticals was derived by Dunlo pet al.l 
d2003l) as M bh = 1.3 x 10 _6 L^ u 3 ^ e , in solar units. We overplot 
this relation as a dashed line in Fig. [6] There is a surprisingly 
good agreement between the low-z relation and the high-z data, 
which at face value might indicate no evolution at all. Of course, 
this cannot be true. 

It is important to recall that applying the above formula ex- 
plicitly assumes a well-defined mass-to-light ratio for QSO host 

[ alaxies, essentially t hat of old elliptical galaxies at z ~ 
Jorgense n et al.fl996l) . Overplotting this relation into a diagram 
showing data obtained at z > 1.5 up to z ~ 3, without account- 
ing at least for passive evolution, is therefore misleading. Only 
spectral information can reveal the amount of galaxy evolution 
required to match the high-redshift results to the local universe. 
Unfortunately, for most of the data shown in Fig. [6] there are 
only single-band measurements availa ble. But since a s ubstantial 
amount of evolution is unavoidable, iPeng et all d2006al) argued 
that even after conservatively correcting for only passive evolu- 
tion with an assumed formation redshift of z = 5, there must be 
a significant offset between the local and high-redshift relations. 
They concluded that at redshifts z ~ 2, the ratio of black hole 
to stellar bulge masses must have been higher than 'today' by at 
least a factor of 4. We explore in the next subsection how our 
colour measurements fit into this picture. 

Using our black hole masses we can also predict Eddington 
luminosities and Eddington ratios e (cf. Table |3J. We find that 
even these high-luminosity QSOs do not exceed Eddington- 
limited accretion, in fact their inferred Eddington ratios are com- 
pletely consistent with those obtained fo r less luminous AGN at 
all redshifts (e.g jKollmeier et al.ll2006l) . 



5.3. Host galaxy colours: Evidence for young stellar 
populations 

For the two QSO hosts detected in both H and K s we could also 
derive rest-frame B — V colours; for the third object we obtained 
a conservative limit on B — V . We find that the colours of our two 



HE 2149-0436 
2 = 3.643 
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wav elength [A] 




4000 5000 6000 7000 
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Fig. 5. Slit spectra of HE 2149-0436, HE 2355-5457 and 
HE 2348-1444. f x is in units of 10~ 16 erg cm' 2 s" 1 A~\ 



objects depend very little on the assumed scale lengths or mor- 
phology. The results are compiled in Table [3] The host galaxies 
are remarkably blue, with a B — V of just below (in the Vega 
system) with an uncertainty of 0.3 mag, including our individ- 
ual uncertainties in the host magnitudes of about 0.2 mag and in 
the K-correction of about 0.05 mag) - thus corresponding to the 
colours of A-type stars. In terms of luminosity-weighted stellar 
ages (the ages of single-burst stellar populations (SSP) with the 
same B — V), we obtain values of ~ 300 Myrs for both objects. 

Finding that QSO host galaxies at z ~ 3 have blue colours 
and young stellar populations is not really surprising, but it nev- 
ertheless poses relevant constraints. Note that a literal interpre- 
tation of the SSP ages leads to a formation redshift of z ~ 3.3, 
considerabl y lower than the z = 5 formation (conservatively) 
assumed bv lPeng et al.1 d2006al) . 
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Table 3. Host galaxy and nuclear properties: Observed H — K s , rest-frame B — V, stellar mass M*, C iv emission line dispersion 
a\,C iv continuum fluxes at 1350A, estimated black hole mass Mbh, Eddington ratio e and the black hole / bulge mass ratio /i. 
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Fig. 6. Relation between black hole masses and host luminosi- 
ties. Filled squares are our data; additional data points are taken 
from the literature. The dotted line shows the zero-redshift rela- 
tion from lDunlop et alj d2003l) . 
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Fig. 7. Mbh-Mbuigo relation reproduced from Haring & Rixl 
(2004) for a sample of local galaxies. The solid line represents 
a constant mass ratio -Mbh/Aftmigc = 0.0014. The three filled 
symbols show our estimated stellar and virial black hole masses 
for the three z ~ 3 QSOs. 



Blue colours of QSO host galaxies hav e recently been estab- 
lished as ubiquitous at lower redshifts. In Jahn ke et al.l d2004a) 
we demonstrated that z < 0.2 quasar hosts with elliptical mor- 
phology have significant blue r colours than their inactive coun- 
terparts. [Sdnchez^etal] (|2004) presented a similar conclusion for 
a sample of QSO hosts at z ~ 0.7 using HST imaging. In a 
different approach, Kauffma nn et al.l (120031) used SDSS data to 
show that Seyfert-2 galaxies are typically hosted by massive 
early-type galaxies with strong evidence for recent star forma- 
tion. Comparing our new rest-frame B — V values with those 
typically found at low redshifts, we find them to be even some- 
what bluer but formally still consistent (given our error bars). 

Corresponding data on high redshift QSO hosts are very 
scarc e. From optical imag ing of several z > 1.8 QSOs with 
HST I Jahnke et al. (2004b) found evidence for substantial young 
stellar populations, with rest-frame UV luminosities compara- 
ble to Lyman Break Galaxies. We do not sample the rest-frame 
UV in our objects (and presumably would be overwhelmed by 
the nuclei anyway), so a direct comparison is not possible. We 
conclude that our finding of blue colours is qualitatively in very 
good agreement with most previous results. 



5.4. Stellar mass-to-light ratios and the M hh -M huigc 
relation 

We can now return to improve our handling of the relation be- 
tween black holes and their host galaxies. While there is no 
unique conversion of a single colour measurement into a stel- 
lar mass-to-light ratio, it is clear that the blu e B — V colours 
found for our objects imply rather low Ml L. iBell et al. (2003) 
demonstrated that for a variety of star formation histories, stellar 
M I L can be quite robustly predicted from B — V, and we used 
their formula 

log w (M/L v ) = -0.628 + (1.305 x (B — V)) . (2) 

where M / Ly is given in solar units. 

Using the B — V colours of Table [3] we obtained M/Ly 
ratios in the range of 0.05-0.2. Given our estimated uncertain- 
ties for the colours the associated uncertainties in M/Ly are 
±0.25 dex 

Combining our M/L ratio estimates with the host galaxy 
luminsities above leads to stellar masses in the range of a few 
10 11 Mq for the two main objects. The formally computed val- 
ues are listed in Table [3] Thus, these galaxies are still very mas- 
sive, but well within the range of the empirically samples stel- 
lar mass function of inactive galaxies as determined recently by 
iFontana et all J2006I) . 
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We then combined the black hole masses derived above with 
our estimated stellar masses. In doing so we assumed our lu- 
minous, massive host galaxies to be dominated by a bulge-like 
component. Our derived ratio fj, between Mbh and Mbuige is 
given in Table [3j the two values and one upper limit are all 
> 0.01. Formally, this is considerably larger than the local ratio 
of ~ 0.0015 obtained for inactive galaxies, although the discrep- 
ancy is only of the order of l-2er given the substantial individual 
error bars. Note that any additional non-bulge component would 
increase the discrepancy. 

We illustrate our results in Fig. |7]by overplotting our high- 
re dshift data poin t s into the local Aibh-Mb u i g c relation obtained 
bv lHaring & Rixl (|2004|) . together with a constant fi = 0.0014 
relation (solid line). Our points are always above the local rela- 
tion, although only marginally. Furthermore, the masses probed 
by our high-luminosity objects are at the extreme end of the 
range probed by the local data, and it is not yet clear that it can 
be extrapolated at constant /i into the high mass regime. 

Measuring Mbh/Mb u i gc for different cosmological epochs, 
is presently a hot topic, and many groups have now en- 
deavoured to estimate fi with a variety of different methods. 
Interestingly, most find evidence for n to increase with look- 
back time. We already mentioned Pen g et all d2006allbl) who ob- 
tained /i 2= 2//i 2 =o *S 4 at z ~ 2, with the inequality depending 
on the unknown stellar mass-t o-light ratios of their galaxies. A 
similar result was obtained by iMcLure et al.l d2006l) . combining 
-Mbuige of radio galaxies with Mbh of radio-loud QSOs, assum- 
ing these to be actually the same class of objects. An extre mely 
large evolution in /j,(z) was claimed by Shi elds et al.l ([2006) who 
used CO emission line widths as surrogates for stellar veloc- 
ity dispersion a nd found that h z= a I 'ij, z =q ~ 50 at z ~ 4. On 
the other hand, lLauer et ail d2007l) recently demonstrated that 
comparing high-redshift quasar measurements with low -redshift 
galaxies can lead to systematically overestimating the amount 
of evolution in /i. Our experimental setup is clearly affected by 
such biases, as the quasars were drawn from a surve y for very lu- 
minous quasars. From Fig. 6 in lLauer et al.l d2007l) we estimate 
that our value of // may be biased by up to 0.5 dex towards high 
values. Taking this bias into account then there is still a system- 
atic offset of our n values to the local one, but it is now well 
below our la errors. 

6. Conclusions 

We have presented rest-frame B and V imaging of QSO host 
galaxies, allowing us for the first time to measure an optical 
B — V colour for at least some of our targets. We found these 
galaxies to show rather, but not extremely blue colours, indica- 
tive of either an enormous burst of star formation a few hundred 
Myr ago, or of ongoing star formation on a much lower but still 
substantial level. Clearly, a single colour cannot differentiate be- 
tween these options. It might be useful to search for CO emis- 
sion in these QSOs which, if found, could be taken as evidence 
for ongoing star formation. 

At least one of our QSOs is embedded in a highly asymmet- 
ric host galaxy structure, indicative of gravitational interaction 
or merging. It is interesting to note that this asymmetric struc- 
ture is more prominent in the H (rest-frame B) than in the K s 
(rest-frame V) band image, which suggests that the blue light 
from young stars comes preferentially from this structure. 

Our observations show that high-luminosity QSOs (Mb ^ 
—28) at z ~ 3 already live in very massive galaxies of several 
times 10 11 M Q . Such massive str uctures are rare at high z, but 
not exceedingly rare: Following lFontana et al.l (|2006), the space 



density of galaxies with stellar masses above 3 x 10 11 M at 
z ~ 3 is of the order of 10~ 5 -10~ 6 Mpc~ 3 , still su bstantially 
highe r than the number density of luminous QSOs dWisotzkil 
2000). Thus, even at the peak of the cosmic QSO activity, only 
a small fraction of galaxies is growing their black holes at near 
Eddington rates. 

By combining stellar bulge masses and virial black hole 
masses for two objects, we could estimate the ratio /i of these 
quantities in individual objects. We found a marginal excess of 
fj,(z — 3) compared to the local relation. Taking into account 
the substantial uncertainties, and considering recent analyses of 
systematic biases, however, our results are still consistent with 
no evolution in the black hole mass - bulge mass ratio. Still we 
believe that the approach as such is promising, and with larger 
samples and more accurate measurements, obtained at several 
redshifts, it has the potential to significantly constrain the degree 
of co-evolution between black holes and hosting bulges. 
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